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Abstract The possibility of Bose-Einstein condensation of excitations has been 
discussed for a long time. The phenomenon of the phase-coherent precession of 
magnetization in superfluid 3 He and the related effects of spin superfluidity are 
based on the true Bose-Einstein condensation of magnons. Several different states 
of coherent precession has been observed in 3 He-B: homogeneously precessing 
domain (HPD); persistent signal formed by Q-balls at very low temperatures; co- 
herent precession with fractional magnetization; and two new modes of the co- 
herent precession in compressed aerogel. In compressed aerogel the coherent pre- 
cession has been also found in 3 He-A. Here we demonstrate that all these cases 
are examples of a Bose-Einstein condensation of magnons, with the magnon in- 
teraction term in the Gross-Pitaevskii equation being provided by different types 
of spin-orbit coupling in the background of the coherent precession. 
PACS numbers: 67.57.Fg, 05.45.Yv, 11.27.+d 



1 Introduction 

Bose-Einstein condensation (BEC) is a phenomenon of formation of collective 
quantum state, in which the macroscopic number of particles is governed by a sin- 
gle wave function. The phenomenon of Bose-Einstein condensate was predicted 
by Einstein in 1925. For a review see, for example Ref. -. The almost perfect 
BEC state was observed in ultra could atomic gases. In Bose liquids, the BEC 
is strongly modified by interactions, but still remains the key mechanism for the 
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formation of a coherent quantum state in Bose systems, which exhibits the phe- 
nomenon of superfluidity characterized by non-dissipative superfiuid mass cur- 
rents discovered first in 4 He by P.L. Kapitza 2 . Superfluidity proved to be a more 
general phenomenon: superfiuid mass current has been found in Fermi liquid 3 He; 
superfluidity of electric charge - superconductivity - is known in metals; quantum 
chromodynamics is characterized by superfluidity of chiral charge; color superflu- 
idity is discussed in quark matter and baryonic superfluidity - in neutron stars; etc. 
Here we discuss the magnon BEC which leads to superfluidity of spin current. 

Strictly speaking, the theory of superfluidity and Bose -Einstein condensation 
is applicable to systems with conserved U(l) charge or particle number. How- 
ever, it can be extended to systems with a weakly violated conservation law. This 
means that it can be applicable to a system of sufficiently long-lived quasiparti- 
cles - discrete quanta of energy that can be treated as condensed matter counter- 
part of elementary particles. In magnetically ordered materials, the corresponding 
propagating excitations are magnons - quanta of spin waves. Under stationary 
conditions the density of thermal magnons is small, but they can be pumped by 
resonance radio-frequency (RF) field (magnetic resonance). One may expect that 
at very low temperatures, the non-equilibrium gas of magnons could live a rela- 
tively long time, sufficient for formation of coherent magnon condensate. 

Recently there appeared a number of articles, where authors claimed the obser- 
vation of BEC of quasiparticles: excitons 3 and magnons^. According to Snoke* to 
claim the observation of BEC one should demonstrate the spontaneous emergence 
of coherence. And, even better, to show the interference between two condensates. 
Since the spontaneous coherence has not been observed directly, the claim of ob- 
servation of BEC in the above article is still under question. 

In superfiuid 3 He-B, the formation of a coherent state of magnons was discov- 
ered about 20 years ago 6 . In pulsed NMR experiments the spontaneous formation 
of domain with fully phase-coherent BEC of magnons has been observed even in 
the presence of inhomogeneous magnetic field. This domain was called the Homo- 
geneously Precessing Domain (HPD). The main feature of HPD is the induction 
decay signal, which rings in many orders of magnitude longer, than prescribed by 
inhomogeneity of magnetic field. This means that spins precess NOT with a local 
Larmor frequency, but precess coherently with a common frequency and phase. 
This BEC can be also created and stabilized by continuous NMR pumping. In 
this case the NMR frequency plays a role of magnon chemical potential, which 
determines the density of magnon condensate. The interference between two con- 
densates has also been demonstrated. It was shown that HPD exhibits all the prop- 
erties of spin superfluidity (see Reviews 7,8 ). The main property is the existence of 
spin supercurrent, which transports the magnetization on a macroscopic distance 
more than 1 cm long. This spin supercurrent flows separately from the mass cur- 
rent, which is essentially different from the spin-polarized 3 He-Ai, where spin is 
transported by the mass current. Also the related phenomena have been observed: 
spin current Josephson effect; phase-slip processes at the critical current; and spin 
current vortex - a topological defect which is the analog of a quantized vortex in 
superfluids and of an Abrikosov vortex in superconductors; etc. 

The spin-orbit coupling, which is responsible for the interaction between magnons, 
is relatively small in 3 He-B. As a result HPD represents almost pure BEC of 
magnons 910 . In typical 3 He experiments, the critical temperature of magnon con- 
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densation is 3 orders of magnitude higher, than the temperature of superfluid tran- 
sition; i.e. magnons undergo the condensation as soon as chemical potential and 
spin-orbit coupling allow for this process. The superfluid 3 He is a very unique 
complex macroscopic quantum system with broken spin, orbital and gauge sym- 
metries, where the spin-orbit coupling due to dipole-dipole interaction between 
the spins of 3 He atoms can be varied experimentally. Under different conditions 
one can observe different types of the BEC of the gas of magnons in 3 He-B— ; 
non-topological solitons called Q-balls in high energy physics^; and also magnon 
BEC in 3 He- A—. 



2 Coherent precession as magnon BEC 

As distinct from the static equilibrium magnetic states with broken symmetry, 
the phase-coherent precession is the dynamical state which experiences the off- 
diagonal long-range order: 

(S+)=S+=S sin Pe im+ia . (1) 

Here S+ is the operator of spin creation; S+ = S x + iS y ; S = (S x ,S y ,S z = Scos/3) 
is the vector of spin density precessing in the applied magnetic field H = Hz; 
P, ft) and a are correspondingly the tipping angle, frequency and the phase of 
precession. In the modes under discussion, the magnitude of the precessing spin 
S equals to an equilibrium value of spin density S = %H jy in the applied field, 
where % is spin susceptibility of 3 He-B, and y the gyromagnetic ratio of the 3 He 
atom (the coherent state with half of magnetization S = (1 /2)%H /yhas been also 
observed in bulk 3 He-B 14 ). Similar to the conventional mass superfluidity which 
also experiences the off-diagonal long-range order, the spin precession in Eq.([T|) 
can be rewritten in terms of the complex scalar order paramete r 7 ^ 10 

(W) = W = y/2Sjhsm | e im+ia , (2) 

If the spin-orbit interaction is small and its contribution to the spectrum of magnons 
is neglected in the main approximation (as it typically occurs in 3 He), then W co- 
incides with the operator of the annihilation of magnons, with the number density 
of magnons being equal to condensate density: 

n M = (^4>) = \V\ 2 = S - 1 ±. (3) 

This implies that the precessing states in superfluid 3 He realize the almost com- 
plete BEC of magnons. The small spin-orbit coupling produces a weak interac- 
tion between magnons and leads to the interaction term in corresponding Gross- 
Pitaevskii equation for the BEC of magnons (further we use units with Ti = 1): 
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Here the role of the chemical potential jJ. = ft) — Ol is played by the shift of the pre- 
cession frequency from the Larmor value (Ol = jH; the latter may slightly depend 
on coordinates if the field gradient is applied. In coherent states, the precession 
frequency ft) is the same throughout the whole sample even in the nonuniform 
field; it is determined by the number of magnons in BEC, Nm = J d 3 rriM, which is 
conserved quantity if the dipole interaction is neglected. In the regime of continu- 
ous NMR, ft) is the frequency of the applied RF field, ft) = (Orf, and the chemical 
potential Li = (Orf — (Ol determines the magnon density. Finally, is the magnon 
mass; and Ed the dipole interaction averaged over the fast precession. The general 
form of Ed(\W\ 2 ) depends on the orientation of the orbital degrees of freedom 
described by the unit vector 1 of the orbital momentum, see Ref.— . 



3 Magnons BEC in bulk 3 He-B and in non-deformed aerogel 

In the coherent precession in bulk 3 He-B and in a non-deformed aerogel—, the 
spin-orbit coupling orients vector 1 along the axis of precession, i.e. 1 || H. In this 
case the interaction term Ed expressed through the condensate order parameter 
has the form different from conventional 4-th order term in dilute gases 9 : 

E D = 0, |ff<^ ; (6) 

>2(m 2 5\ 2 lu , l2 . 5. 



E D = -xnt^--) ,\Vr>-S. (7) 

Here Ql <C (Ol is the Leggett frequency which characterizes the dipole interac- 
tion. If the chemical potential ji is negative, i.e. ft) is less than (Ol, the minimum 
of the Ginzburg-Landau (GL) energy Ed(\ 1 I , \ 2 ) — jx\W\ 2 corresponds to f = 0, 
i.e. to the static state with non-precessing equilibrium magnetization (j3 = 0). For 
jU > (i.e. for positive frequency shift), the profile of the (normalized) energy 
density is shown in Fig. [TJ for several values of jj. given in dimensionless units 
A = (<°rf — ~ I 1 /(^l/^ (0 l)- The minimum of the GL energy corre- 

sponds to \*F\ 2 /S = (5/4) + (15/32)jtt. The consequence of the peculiar profile 
of the interaction term is that as distinct from the dilute gases the formation of the 
magnon BEC starts with the finite magnitude \ W\ 2 = (5/4)5. This means that the 
coherent precession starts with a tipping angle equal to the magic Leggett angle, 
P = 104°, and then the tipping angle increases with increasing frequency shift. 
This coherent state called the HPD persists indefinitely, if one applies a small RF 
field to compensate the losses of magnons caused by small spin-orbit interaction. 

In conventional magnetic systems, magnetization precesses in the local field 
with the local frequency shift and thus experiences dephasing in the inhomoge- 
neous field. In the case of magnon BEC, the rigidity of the order parameter (the 
gradient term in Eq.[5]l plays an important role. The spatial dephasing leads to the 
gradient of chemical potential. This in turn excites the spin supercurrents, which 
finally equilibrate the chemical potential. In the steady state of magnon BEC the 
gradient of a local field is compensated by small gradient of magnon density \ W\ 2 
in such a way, that the precession frequency and its phase remain homogeneous 
throughout the whole sample. 
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Fig. 1 (Color online) The Ginzburg-Landau energy Eo(\ t I'\ 2 ) — /x|f| 2 in bulk 3 He-B as a func- 
tion of magnon density (tipping angle of precession) for different values of the normalized chem- 
ical potential fi = (£i)| F — a>l)/i2l, Energy minima correspond to magnon BEC, i.e. coherent 
HPD states precessing with frequency shift equal to the chemical potential 
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Fig. 2 (Color online) The typical signal of induction decay from the BEC of magnons; (a) stro- 
boscopic record of the signal; (b) amplitude of signal; (c) frequency of the signal. Frequencies 
469.95 and 469.4 kHz correspond to Larmor frequency at the top and the bottom of the cell. 
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Fig. 3 (Color online) The Ginzburg-Landau energy E D ( | f | 2 ) - ii | f | 2 as a function of magnon 
density for the same frequency to, but at different Larmor frequency (Ol at the top and the bottom 
of the cell (this corresponds to different normalized chemical potential fi = (tt) 2 — mj)/Ql) 



In a pulsed NMR experiment, the magnetization is deflected by a strong RF 
pulse. The typical induction signal after the pulse in the cell with a large gradient 
of magnetic field along the axis of the cell is shown in Fig. [2] Due to the field gra- 
dient the induction signal should dephase and disappear in about 10 ms. Instead, 
after a transient process of about 2 ms, the induction signal acquires an amplitude 
corresponding to a 100% coherent precession of the deflected magnetization with 
the spontaneously emerging phase. This coherent state lives 500 times longer than 
the dephasing time caused by inhomogeneity. 

What happens with the magnon system during the transient period? The GL 
energy Ed(\W\ 2 ) — /i|f| 2 in the presence of the gradient of magnetic field is 
shown in Fig. [3] As soon, as at the higher field end of the cell the deflection of 
magnetization becomes smaller than 104°, the spin-orbit interaction cannot any- 
more compensate the gradient of magnetic field. The cell splits into two domains, 
in one of them the magnetization is stationary, while in other one magnons con- 
dense with the density close to |f| 2 = (5/4)5, i.e. the magnetization precesses 
with j5 slightly above 104°. In the subsequent process of relaxation caused by the 
non-conservation of magnon number, the volume of the BEC condensate (HPD) 
decreases. During the relaxation, the BEC does not loose the phase coherence, but 
its chemical potential (precession frequency) changes. The frequency of HPD cor- 
responds to Larmor frequency at the boundary of the domain (Fig. [3) and slowly 
changes with relaxation as the boundary moves down (see Fig.|2]l. The amplitude 
of the signal exactly corresponds to the record of the frequency. 

Of course, the HPD has been observed, studied and explained on the basis of 
theory of spin superfluidity and non-linear NMR long time ago 6 . However, the 
consideration of this phenomenon in terms of the magnon BEC not only demon- 
strates the real system with the BEC of excitations, but also allows us to simplify 
the problem and to study and search for the other types of the magnon BEC in 3 He, 
such as Q-balls 12 ; HPD2 state found in a deformed aerogel 1 1 ; coherent precession 
in 3 He- A also found in the deformed aerogel 13 ; etc. 
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Fig. 4 (Color online) The normalized GL energy Sodfl 2 ) - At|f I 2 in superfluid 3 He-A in the 
deformed aerogel, which orients the orbital momentum 1 parallel to the magnetic field. The GL 
energy is given as a function of magnon density (tipping angle of precession) for different values 
of the normalized chemical potential tx = (to 2 — a>l) The minima of GL energy correspond 
to the stable magnon BEC 



4 Magnons BEC in superfluid 3 He in aerogel 

In superfluid 3 He- A, the spin-orbit interaction depends on the polar angle of 
the orbital vector 1 in the following way: 



Ed 



(8) 



While the sign of the quadratic term in Eq. © is not important because this term 
can be compensated by the shift of the chemical potential fx in Eq. [5] the sign 
of the quartic term is crucial for the stability of the BEC. In the static state (i.e. 
f — 0) of a bulk 3 He-A, the dipole energy is minimized when 1 is perpendicular to 
the magnetic field, sin = 1 . In NMR, this orientation of 1 remains frozen, and the 
quartic term remains negative. This attractive interaction between magnons desta- 
bilizes the BEC, which means that homogeneous precession of magnetization in 
3 He- A is unstable, as was predicted by Fomin 17 and observed experimentally—. 

Recently, we found the method, how to reorient the orbital momentum 1. For 
that we immersed 3 He inside the squeezed aerogel. Even a small deformation of 
a few percent fixes the orbital momentum along the deformation. 19 . By applying 
the deformation along the magnetic field we fixed Pi = 0, which made positive the 
quartic term in Eq. ([8). As a results magnons can condense into the stable BEC 
state, as we observed experimentally—. For /3l = 0, the GL energy Eo(\W\ 2 ) — 
/I |f| 2 (Fig. |4} is similar to the traditional GL energy for the conventional BEC in 
atomic gases. This means that, as distinct from the HPD in 3 He-B, the minimum of 
the GL energy is realized even for small magnon densities, and in the experiments 
made in the Institute of Solid State Physics in Kashiva, Japan, we observed the 
formation of BEC in 3 He-A without the domain formation. 
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Fig. 5 (Color online) The angle of deflection of orbital momentum vector 1 (left scale) and 
magnon density (right scale) as a function of position in the one dimensional cell for different 
total number of magnons Q. The NMR radiation frequency is shown for each curve. 

In the 2007 experiments in Neel Institute, Grenoble we applied the aerogel 
deformation to reorient the orbital momentum 1 in 3 He-B perpendicular to mag- 
netic field. For this case the shape of the GL energy allows for stable and unstable 
magnon BEC depending on the tipping angle f}. Both BEC states were observed 
in pulsed NMR experiments in a good agreement with the GL energy—. 



5 Magnon condensation into Q-ball 

The non-topological soliton solutions, called Q-balls, have been proposed by Cole- 
man— as a mechanism of formation of composite particles in high energy physics. 
Q-balls appear in theories containing a complex scalar field with suitable self- 
interaction. They are stabilized due to the conservation of the particle number, or 
in general the charge Q. The Q-ball is a rather general physical object (see- 2 ^), 
which in principle can be formed in different condensed matter systems. In partic- 
ular, Q-balls were suggested in the atomic Bose-Einstein condensates^ 2 -. 

The Q-balls were found in NMR experiments in superfiuid 3 He-B as a small- 
amplitude, but extremely long lived induction decay signal 23 . The role of the Q- 
charge is played by the total number of magnons Q = J cPr um—- In all previous 
examples, magnons were almost homogeneously distributed inside the BEC state. 
In Q-balls, the condensation of magnons occurs in a local trap, produced by the 
texture of orbital momentum 1. In the homogeneous state with 1 || H the spin-orbit 
interaction is zero for small W (see Eq. (O), but in the 1-texture Ed ^ 0: 



Ed = X&l 



2 r 4sin 2 (fe/2) 
L 5S 



2 



sin 4 (fe/2) 4 

S 2 \ r \ 



(9) 



where Pi is the polar angle of the vector 1. It contains the negative quartic term, 
which describes the attractive interaction between magnons. This form does not 
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support the Bose condensation in bulk 3 He, but supports formation of magnon 
liquid droplets. Such a droplet may form inside of the existing orbital texture trap. 
Also, the self-localized orbital trap can be spontaneously created by magnons, in a 
full analogy with Q-balls. The magnon liquid in the droplet and inside the Q-ball 
remains in the coherent precessing state. 

The droplet or Q-ball is formed in places of the sample where the potential 

5m L 2 

produced by the 1-texture has a minimum. The precise form of a Q-ball depends on 
the particular texture and on the position in the container. Typically it is formed in 
the so-called flared-out texture realized in cylindrical cells. Far from the horizontal 
walls and close to the axis of the cell the angle fii linearly depends on the distance 
r from the axis: Pi(r) « fcr (see review—). The potential d 1 Ob for magnons is 
U(r) oc K 2 r 2 , similar to the harmonic trap used for confinement of dilute Bose 
gases. The peculiarity of the attractive quartic term in the GL functional in Eq. (O 
is that it is not a constant but is oc jc 4 /- 4 . This stabilizes the condensate in the trap: 
at fixed number Q of magnons the GL energy has minimum at some Q-dependent 
size of the magnon droplet, r(Q). This is the main feature of Q-balls. Fig. [5] 
shows numerical simulations of the spatial distribution of magnon density inside 
the 1 texture for different Q in one dimensional case. Note another feature of Q- 
balls: with increasing number Q of the trapped magnons the potential well itself 
becomes modified by magnons and finally it is overwhelmingly determined by Q. 

Simultanious formation of many different Q-balls has been observed in CW 
NMR 25 . They are generated by spin waves in the broad range of frequencies. 
They can be created in particular by topological defects at the walls of the cell. 
As distinct from CW NMR, where the Q-balls are generated starting continuously 
from Q = 0, in pulsed NMR a Q-ball is formed after a large Q is pumped into the 
cell. In this case the 3D Q-ball is often formed on the axis of the flared-out texture, 
away from the horizontal walls 26 . Finally the lowest energy Q-ball collects all the 
magnons and radiates as a single harmonic oscillator with a slowly increasing 
frequency, due to the modification of the trap potential in the process of decay of 
the Q charge, as shown in Fig. [5] The relaxation of lowest energy Q-balls can be 
compensated by additional RF pumping of magnons at any frequency above the 
frequency of Q-ball, as was found in 27 and explained in—. 



6 Conclusion 

The variety of observed BEC states of coherent precession in superfluid 3 He is 
provided by specific interaction terms in the GL energy, which in general are dif- 
ferent from the conventional 4-th order term in the traditional BEC. They come 
from the spin-orbit interaction and are governed by the orbital vector 1. This al- 
lows us to regulate the profile of the GL energy: by applied counterflow 29 ; by 
deformation of aerogel 1311 ; by exciting the precession at one half of equilibrium 
magnetization 14 ; by using flared-out 1-texture 12 and topological defects; etc. 
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